Geometric measure of quantum discord and the geometry of a class of two-qubit states 



O 
(N 

o 

Q 



Wei Song/ Long-Bao Yu,^ Ping Dong/ Da-Chuang Li/'^ Ming Yang,'' and Zhuo-Liang Cao^ 

^School of Electronic and Information Engineering, Hefei Normal University, Hefei 230061, China 
^Hefei National Laboratory for Physical Sciences at the Microscale and Department of Modern Physics, 
University of Science and Technology of China, Hefei 230026, China 
'^School of Physics and Material Science, Anhui University, Hefei 230039, China 

(Dated: December 22, 2011) 

We investigate the geometric picture of the level surface of quantum entanglement and geometric 
measure of quantum discord(GMQD) of a class of X states, respectively. This pictorial approach 
provides us a direct understanding of the structure of entanglement and GMQD. The dynamic 
evolution of GMQD under two typical kinds of quantum decoherence channels is also investigated. 
It is shown that there exist a class of initial states for which the GMQD is not destroyed by 
decoherence in a finite time interval. Furthermore, we establish a factorization law between the 
initial and final GMQD which allows us to infer the evolution of discord under the influences of the 
environment. 
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I. INTRODUCTION 

Entanglement is regarded as an ingredient resource for 
performing almost all quantum information processing 
tasks[l|. The situation started to change until a compu- 
tational model was presented named deterministic quan- 
tum computation with one qubit(DQCl)|i|. Quantum 
discord was considered to be the figure of merit for this 
model of quantum computation. Ever since, quantum 
discord has attracted much attention (3l427| . Quantum 
discord is a measure of nonclassical correlations between 
two subsystems of a quantum system. It quantifies how 
much a system can be disturbed when people observe it 
to gather information. Such quantum correlations may 
be present in separable states and have a non vanish- 
ing value for almost all quantum state [l^. On the other 
hand, the quantum discord can be used to indicate the 
quantum phase transitions better than entanglement in 
certain physical systems at finite temperature (28l . [29j . 
In particular, two operational interpretations of quantum 
discord have been proposed, one in thermodynamics (30| 
and the other from the information theoretic perspective 
through the state merging protocol (sil . [3^ . These re- 
sults establish the status of quantum discord as another 
important resource for quantum informational processing 
tasks besides entanglement. 

Originally, the first definition of quantum discord is 
due to Ollivier and Zurek Qand, independently, by Hen- 
derson and Vedral However, it is notoriously difficult 
to compute it because of the minimization taken over 
all possible POVM, or von Neumann measurements. At 
present, there are only a few analytical results includ- 
ing the Bell-diagonal states [Bj, rank-2 states [Tol. [33|and 
Gaussian state [S^ . [ssj . In addition, a simple algorithm 
to evaluate the quantum discord for two-qubit X-states 
is proposed by Ali et al. (36| with minimization taken 
over only a few cases. Unfortunately, their algorithm 
is valid only for a family of X-states d, [s^l- Recently 
Shi et al. [3^ present an efficient method to solve this 



problem. For the general two-qubit states, the evalua- 
tion of quantum discord remains a nontrivial task and 
only some lower and upper bounds are available (39j . 
In order to avoid the difficulties in minimization proce- 
dures a geometric view of quantum discord was intro- 
duced. Generally, there are two version of geometric 
measure of quantum discord(GMQD). In the first ver- 
sion the concept of relative entropy is used as a distance 
measure of correlations Q. The second version is de- 
fined by the Hilbert-Schmidt norm measure 0. The 
relative-entropy-based discord have the drawback that 
their analytical expressions are known only for certain 
limited classes of states. Below we only consider the sec- 
ond version of geometric measure of quantum discord. 
Formerly, this geometric measure of quantum discord is 
defined by D\ = min t|p — xll^i where fig denotes the 

set of zero- discord states and \\X - Y\\ = Tr {X - Y) 
is the square norm in the Hilbert-Schmidt space. The 
subscript A denotes that the measurement is taken on 
the system A. An arbitrary two-qubit state can be writ- 
ten in Bloch representation: 
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I ® ffi) + RijCTi (8) (Tj I (1) 



where Xi = Trp{ai ® I) ,yi = Trp{I ® Oi) are compo- 
nents of the local Bloch vectors, ai,i € {1,2,3} are the 
three Pauli matrices, and Rij ~ Trpiai cTj) are com- 
ponents of the correlation tensor. For two-qubit case, 
the zero-discord state is of the form X = Pi\4'i) ("01 1 ® 
Pi+P2\ip2) (i/'2| ® P2, where {\^pl) ,\ip2)} T-S a. single-qubit 
orthonormal basis. Then a analytic expression of the 
GMQD is given by 0: 
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where x = {xi , 12 , X3) and fcmax is the largest eigenvalue 
of matrix K = xx^ + RR^ . 

In this paper we shall investigate the level surfaces of 
entanglement and GMQD for a class of two-qubit states 
using the geometric picture presented in [4(]j . It is well 
known that the set of Bell-diagonal states for two qubits 
can be depicted as a tetrahedron in three dimensions in 
Bloch representation [4l|. Analogous to entanglement, 
Lang and Caves have depicted the level surfaces of quan- 
tum discord for Bell-diagonal states. They found that the 
picture and the quantum discord are very different. This 
pictorial approach provides a complete and direct under- 
standing of the behavior of von-Neumann measurement- 
based quantum discord. It has been shown that the quan- 
tum discord and GMQD do not necessarily imply the 
same ordering of two-qubit X-states |42| - |44| . More re- 
cently, Girolami and Adesso [i^ and independently Ba- 
tle et al. (46| provided a numerical evidences, from which 
one can infer that there exist other states violating the 
states ordering with quantum discords. It is thus worth 
investigating the GMQD from the geometric picture. In 
this sense, our research provides a direct understanding of 
the structure of GMQD at first time. We consider a class 
of X states that the Bloch vectors are z directional, which 
including Bell-diagonal states as a special case. We show 
that the level surface of GMQD is very different with the 
von-Neumann-nieasurement-based quantum discord. On 
the other hand, the dynamics of quantum discord has at- 
tracted much attention due to the inevitable interaction 
with environment. The dynamic behavior of GMQD is 
also investigated under two typical kinds of decoherence 
channels. We find a class of states for which the GMQD is 
not destroyed by decoherence in a finite time interval. In- 
terestingly, we also obtain a factorization law for GMQD 
which allows us to infer the evolution of discord under 
the influences of the environment, e.g. the phase damp- 
ing channel(PDC), and the depolarizing channel(DPC). 

The paper is organized as follows. In Sec. II, we present 
the geometric picture of the level surface of quantum en- 
tanglement and GMQD of a class of X states with the 
Bloch vectors are z directional, respectively. It is shown 
that the surface of constant GMQD varies with the lo- 
cal Bloch vectors. In Sec. Ill, we investigate the dynamic 
evolution of GMQD under quantum decoherence chan- 
nels, and obtain analytic results for two typical kinds of 
quantum decoherence channels. A summary is given in 
Sec.IV. 



II. GEOMETRICAL PICTURE OF 
ENTANGLEMENT AND GMQD 

For analytical simplicity, we consider the following two- 
qubit X states: 
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where we choose the Bloch vectors are z directional with 
r = (0,0,r),s = (0,0, s). The GMQD can be calcu- 
lated explicitly for this state, thus allowing us to get 
analytic results. If r = s = 0, p reduces to the two- 
qubit Bell-diagonal states. Horodecki have shown that 
Bell-diagonal states belongs to a tetrahedron with ver- 
tices (l,-l,l),(-l,l,l),(l,l,-l),and (-1,-1,-1) in 
the Bloch representation. From the positively of the 
eigenvalues of p in Eq.(3), we have 
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For fixed parameters r and s, the above inequalities be- 
comes a three-parameter set, whose geometry can be de- 
picted in the three dimensional correlation state space. 
In Fig.l we plot the physical region with different r 
and s, respectively. Fig.l shows that physical regions 
of the state p shrink with larger r and s. We plot 
in Fig. 2 the level surfaces of constant concurrence with 
fixed r and s for three cases. Here, we choose con- 
currence to measure entanglement which is defined as 
C = max {0, Al — A2 — A3 — A4}, where the A,; are, in 
decreasing order, the square roots of the eigenvalues of 
the matrix pay (g ayp*ay ay where p* is the complex 
conjugate of p. As shown in Fig. 2, the level surfaces of 
constant concurrence for the state p defined in Eq.(3) 
consist of four discrete pieces, and the areas decrease 
with larger concurrence. A extremal case is the four ver- 
tices (1,-1,1), (-1,1,1), (1,1,-1), and (-1,-1,-1) of 
the tetrahedron correspond to the four Bell states with 
maximal concurrence. Finally, we investigate the GMQD 
from the geometric picture. For the state p, the GMQD 
can be calculated in the method presented in Ref. |47| . 
By introducing a matric TZ defined by 
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and 3x4 matric TZ' through deleting the first row of TZ. 
Then the GMQD is given by 



K if) 



max A2 

k 



(6) 



where A^ is the singular values of TZ' . For the two-qubit 
state p in Eq.(3), we obtain 



{4 -^4 + 4 + r'- Max [4,4,4 + r^))(J) 
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The geometric picture is depicted in terms of the con- 
stant GMQD in Fig. 3. From these plots one can see that 
the shape of the constant GMQD is quite different from 
von-Neumann-measurement-based quantum discord. It 
also shows different shape for different local Bloch vec- 
tors. The constant surfaces are cut of by the physical 
region of state p. For small discord the surface is contin- 
uous, and it becomes discrete pieces for larger discord. 
At the four vertices (1, -1, 1), (-1, 1, 1), (1, 1, -l),and 
(— 1,— 1,— 1) of the tetrahedron the GMQD reach its 
maximal value. Furthermore, one can see that GMQD 
is neither concave nor convex as shown in Fig. 3. 



III. DYNAMICS OF GMQD UNDER LOCAL 
DECOHERENCE CHANNELS 



In this section we consider the state is affected by 
the action of two independent channels and calculate 
the GMQD analytically. The dynamics of quantum dis- 
cord has been investigated in both Markovian and non- 
Markovian environments and has been demonstrated ex- 
perimentally [48l - l56l |. It has been shown that the behav- 
iors of von-Neumann-measurement-based quantum dis- 
cord and GMQD may be different, it is thus desirable to 
consider the evolution of GMQD under different decoher- 
ence channels. Here, we consider two typical kinds of de- 
coherence channels: the phase damping channel(PDC); 
and the depolarizing channel(DPC). To calculate the dy- 
namics of GMQD, we turn to the Heisenberg picture to 
describe the quantum channels. In order to obtain the 
analytic expressions of GMQD of the state subject to 
local decoherence channels, we need to calculate the ex- 
pection matrix TZ. In the Heisenberg picture [U [13 , the 
expectation matrix TZ is given by 



(8) 



where TZq = 
initial state, 



Tr {oi ® UjPQ), i G {0,1,2,3}, po is the 



and M 



A(B) 



is the transmission matrix of 



each local channel. For simplicity, we choose the local 
channels to be identical. In this case, the transmission 
matrixs can be written as 
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For state (3), 7^o = {{1, 0, 0, s} , {0, ci, 0, 0} , {0, 0, C2, 0} 
{r, 0, 0, C3}}. According to the above formula, we have 
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For some Bell-diagonal states, it has been shown 
that von-Neumann-measurement-based quantum discord 
is not destroyed by decoherence for some finite time in- 
terval [5l|. A natural question arises: Whether such 
phenomena exists for GMQD? We consider the state 
in Eq.(3) undergoes two identical PDCs. In this case, 
p = 1 — exp(— 7t), where 7 is the phase damping rate. 



For ci = 0,C2 > 



Suppose pi satisfy the equa- 



tion (1 — pi) C2 = + Cg. If p < pi, from Eq.(ll) we 
have Dpjj^j = ^ (r^ 4- c§) which is independent of time. 
Therefore, we conclude that for a finite time interval 
the GMQD does not decay despite the presence of lo- 
cal phase damping noises. It is directly to see that such 
phenomenon also exists for the case C2 = 0, > + c^. 
These results show that GMQD remain intact under the 
action of some special kinds of quantum channels. It 
should be noted similar phenomenon have also been no- 
ticed by Karpat et al j59| for the qubit-qutrit systems. In 
the geometric picture, this behavior corresponds to the 
state evolves along a straight line in the constant GMQD 
tube until it enconters another constant GMQD tube. 

Hitherto, the time evolution of GMQD under PDC or 
DPC can be described by Eq.(ll) and Eq.(12). Next, we 
want to derive a more general result on GMQD relating 
the initial and final state of GMQD. Inspired by the fa- 
mous factorization Law for entanglement decay derived 
by Konrad et al. we find a analogous factorization 
law between the initial and final GMQD of the class of 
two-qubit states defined in Eq.(3) subject to two differ- 
ent local decoherence channels. We generalize our result 
as the following theorem. 

Theorem. Consider the class of X states defined in 
Eq.(3), with each qubit being subject to the local deco- 
herence channels, i.e. the phase damping channel(PDC) 
or the depolarizing channel(DPC). The time evolution of 
GMQD satisfies 



® $2) p{t)] ^ 2D0 [($1 ® $2) |A> {fi.W b(0)](13) 



where the local decoherence channels is represented by 
$1 and $2, po is the initial state and |/3i),i £ {1,2,3,4} 
denotes one of the four Bell states. 

Proof. First we consider the state p(0) is effected by the 
action of two identical local PDCs, the time evolution of 
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GMQD is given by Eq.(ll). For convenience we divide 
the proof into three cases. 

Case 1. + c| > {cf , c|}. In this case, it is directly to 
show that the inequahty becomes equahty. 
Case 2. cf ^ c| > + c§. First, using Eq.(7) we 
have (p(0)) ^ j {^2 + ^'^ + cf) . Suppose po satisfy the 
equation (1 - po f c\ = r"^ + c\. If p ^ po, the GMQD of 
the state p is given by 
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If p > Poi then 
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(15) 

Case 3. ^ + c| ^ c^. For p ^ po or p > po, the 
proof is similar to case 2. 

By simply exchanging ci and C2 we can verify the above 
relations also hold for the cases c| ^ cf ^ r^+Cg and ^ 
r^ + c| ^ c^. For the case of two identical local DPCs, one 
can prove the above results by the same way as proving 
PDCs. In the following, we consider the first qubit is 
subject to the PDC and the second qubit is subject to 
the DPC. The expectation matrix TZ can be calculated 
according to the formula TZ = MaTZ-qM]^, where Af^(B) 
is the transformation matrix of PDC (DPC). Thus, the 
GMQD is given by 



{p{t)) = \[{l-p)'cl + {l^p)'cl+r^ + {l^p)^cl 



\{i-p)'cUi-p)\iy + {i-p)'ci\ 



(Pit)) 
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= 2D3 [{%pDc ® %DPc) (lA) {m [p (0)) (19) 

Case 3. c\ + c\ c^. In this case, the proof is sim- 
ilar to case 2. 

One can also directly verify that the above relation 
holds for the other cases. ■ 

The theorem above provides us a method to compute 
the lower bound of the time evolution of a class of X 
states under two typical kinds of decoherence channels, 
without resort to the time evolution of the underlying 
quantum state itself. This inequality also holds for the 
one-sided PDC or DPC, which is summarized as the fol- 
lowing Corollary: 

Corollary 1. For the class of X states defined in Eq.(3), 
with one qubit being subject to PDC or DPC, we have 



D<> {p (t)) ^ 2DS [($, ® /) (1/3,) (Al)] (p (0)) (20) 

Moreover, for r = s = 0, i.e. Bell diagonal states. Us- 
ing Eq.(12), one can directly calculate the above inequal- 
ity becomes equality which we summarize as follows: 
Corollary 2. For arbitrary Bell-diagonal states subject 
to two DPCs. the evolution of GMQD is given by 
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where we have assumed the parameter p are same in the 
two decoherence channels. Then it suffices to consider 
three separate cases. 
Case 1. +cl^ {c?,c^}. Then 



where the two DPCs may be different. By far we have 
considered the time evolution of GMQD in the presence 
of two typical kinds of local decoherence noise. Another 
important decoherence noise is the amplitude damping 
channel(ADC). There are evidences to show that the 
above relation also hold for two identical ADCs, i.e. 



{p{t))^^[{l~p)\l + {l-p)^c: 
= 2DS [{$pnc ® %DPc) im im Ds (p(0)) (17) 
Case 2. ^ ^ + c-^. Suppose pg satisfies 
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(1 — Po)'^ cl ~ + (1 — Po)^ c§. If p ^ Pg, we have 
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(22) 



In Fig. 4, we plot the evolution of GMQD under two 
identical ADCs and its lower bound in Eq.(22) when 
(a)ci = 0.1, C2 = 0.1, cg = 0.2, r = s = 0.3; (b)ci = 
0.2, C2 = 0.05, C3 = 0.3, r = 0.4, s = 0.1. 
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IV. DISCUSSIONS AND CONCLUSIONS 

In this work, we have investigated the level surfaces of 
GMQD for a class of two-qubit X states from the geo- 
metric picture. First, we plot the physical region for a 
class of two-qubit X states with fixed local Bloch vectors. 
It is shown that physical regions of the state have differ- 
ent geometry with the Bell-diagonal states and shrink 
with larger Bloch vectors. Second, the geometric pic- 
ture is depicted in terms of the constant concurrence and 
GMQD, respectively. We find that the shape of the sur- 
faces has close relationship with the value of GMQD and 
local Bloch vectors. Finally, we also investigate the dy- 
namics of GMQD under two typical kinds of decoherence 
channels and obtain analytic results of the evolution of 
GMQD. It is shown that there exist a class of initial states 
for which the GMQD is not destroyed by decoherence in 
a finite time interval. Moreover, a direct factorization 
relationship between the initial and final GMQD subject 
to two typical kinds of decoherence channels is derived. 
This factorization law allows us to infer the evolution 
of entanglement under the influences of the environment 
without resort to the time evolution of the initial quan- 



tum state itself. An open question is whether this law 
holds under general local decoherence channels. Our re- 
sults imply that fiurther study on the dynamics of GMQD 
is required. 
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FIG. 1: (Color online). The geometry of the set 
of valid states with different r and a, repectively. 
(a)r=s=0.3,(b)r=s=0.5,(c)r=0.4,s=0.1 
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FIG. 2: (Color online). Sur- 

faces of constant concurrence. (a)r=s 
=0.3,C(p)=0.03,(b)r=s=0.5,C(p)=0.35,(c)r=0.4,s=0.1,C(p)=0.03 



9 



1.0 




1.0 



FIG. 3: (Color online). Surfaces of constant GMQD. (a)r=s 
=0.3,L»(p)=0.03,(b)r=s=0.5,D(p)=0.35,(c)r=0.4,s=0.1,D(p)=0.08 
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FIG. 4: (Color online) .Plots of the dynamics of GMQD under 
two identical ADCs(blue line) and its lower bound defined in 
Eq.(22)(red line).(a)ci = 0.1, C2 = 0.1, C3 = 0.2, r = s = 0.3; 
(b)ci = 0.2, C2 = 0.05, C3 = 0.3, r = 0.4, s = 0.1. 



